Dynamics of a passive sliding particle on a randomly fluctuating surface 
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We study the motion of a particle sliding under the action of an external field on a stochastically 
fluctuating one- dimensional Edwards- Wilkinson surface. Numerical simulations using the single- 
step model shows that the mean-square displacement of the sliding particle shows distinct dynamic 
scaling behavior, depending on whether the surface fluctuates faster or slower than the motion of 
the particle. When the surface fluctuations occur on a time scale much smaller than the particle 
motion, we find that the characteristic length scale shows anomalous diffusion with £^{t) ~ t^'^ , where 
(j> « 0.67 from numerical data. On the other hand, when the particle moves faster than the surface, 
its dynamics is controlled by the surface fluctuations and ^(f) ~ . A self-consistent approximation 
predicts that the anomalous diffusion exponent is (/> = § , in good agreement with simulation results. 
We also discuss the possibility of a slow cross-over towards asymptotic diffusive behavior. The 
probability distribution of the displacement has a Gaussian form in both the cases. 



I. INTRODUCTION 



The advection of a passive scalar field (such as temper- 
ature) by a turbulent fluid is a well-known problem [1], 
and is an example of a coupled semi-autonomous system, 
where one of the fields evolves on it own, but affects the 
dynamics of the other. Many such systems have been 
studied in recent times, such as the dynamics of a par- 
ticle passively sliding on a randomly fluctuating surface 
[2-4] , phase separation on a rough substrate [4] and clus- 
tering of particles on a stochastically fluctuating surface 
under the influence of gravity [5]. These studies have 
shown that such systems possess a number of non-trivial 
features. For example, hard-core particles cluster into a 
fluctuation-dominated phase separated state [5] and non- 
interacting particles have non-trivial density correlations 
[2] in the steady state. It is thus highly desirable to aim 
for an understanding of the relation between the time 
evolution of the stochastic fluctuations of the underlying 
field and the dynamics of the passive scalar, and this is 
our primary motivation. 

In this paper, we study the dynamics of a passive 
sliding particle moving on a randomly fluctuating sur- 
face. For concreteness, we choose the surface to be the 
one-dimensional Edwards- Wilkinson surface with uncor- 
related noise [6] , whose steady state and dynamical prop- 
erties are well-known. As we shall see, even this simple 
example has unexpectedly rich properties. 

For clarity of presentation, we shall state the prob- 
lem here. Let us consider a one-dimensional fluctuating 
Edwards- Wilkinson surface, whose equation of motion 
has the form 



dh „5 



T] ; {r]{x, t)Tj{x', t')) = 2D5{x - x')6{t - t') 

(1) 

The Langevin equation for the motion of a particle 
sliding on the surface has the form 



dX{t) 
dt 



-x{t) 



(2) 



where X{t) is the position of the particle on the surface 
and r is a friction coefficient. In principle, one could in- 
clude a random noise term for the motion of the particle 
also, but since its only effect is to add a diffusive term 
to the mean square displacement, we shall drop it from 
further calculations. 

The time evolution of the mean-square displacement 
of the sliding particle is given by the equation, 

ft 



dt 



{X\t)) = 2T^ / dt'^t,t') 



(3) 



where 



= {d,h[xit),t]d.,h[xit'),t'] 



(4) 

is the 'effective' noise correlator for the motion of the 
sliding particle, and is equal to the slope correlator of the 
fiuctuating surface evaluated along the particle trajecto- 
ries. The angular brackets on both side represents aver- 
age over different realizations of the noise r]{x, t). Clearly, 
this is a highly non-trivial quantity to compute, since it 
is not obvious a priori how the slope correlator would be- 
have when averaged over such a set of highly restricted 
paths. 

We outline our main results here. Numerical simula- 
tions show that when the particle moves much faster than 
the surface, its dynamics is controlled by the surface fluc- 
tuations, and the mean square displacement behaves as 
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{X'^{t)) ~ vt. An argument using the concept of over- 
turning valleys confirms this result. On the other hand, 
when the surface fluctuation is fast, the particle moves in 
a very dynamic landscape, and the mean square displace- 
ment shows anomalous diffusion with {X'^{t)) ~ t^"^ over 
the time scales of simulations. A self-consistent approxi- 
mation which assumes no correlation between a particle 
trajectory and the underlying surface height configura- 
tion, predicts that the anomalous diffusion exponent is 

= I , in reasonably good agreement with numerical re- 
sults. We also discuss the possibility that the observed 
anomalous diffusion might be simply transient behavior, 
and give reasons for the same. 

The rest of this paper is arranged in the following way. 
In the next section, we present the self-consistent approx- 
imation to evaluate the correlator in Eq.4, and discuss its 
predictions and limitations. In Sec. Ill, we present the 
results of numerical simulations of the problem using the 
single-step algorithm. In Sec. IV, we summarize our 
findings and discuss the outstanding questions. 



II. THE SELF-CONSISTENT APPROXIMATION 

Let us define z = X{t)—X{t'), which, like X{t) itself, is 
a random variable. Let us now sort the set of all available 
surface configurations using z, i.e., for fixed t and t' , let 
Sz denote the set of surface configurations where a sliding 
particle is displaced by a distance z ± Az between times 
t and t'. The fraction of such configurations is simply 
P{t,t' , z)Az, where P{t,t\z) is the (unknown) proba- 
bility distribution function for particle displacements be- 
tween times t and t'. It then follows that 

/oo 
dzP{t, t', z){d,h[X{t),t]d,h[X{t) + z, t'])s^ 
-oo 

(5) 

Note that the average in Eq.5 is a restricted avc;rage, 
i.e., it is evaluated over surface configurations in which 
the sliding particle is moved by a distance z between t 
and t' . 

Up to now, our treatment has been exact. In order 
to make further progress, we make an approximation by 
assuming that the restricted slope correlator in Eq.5 is 
simply the standard slope correlator of the surface, evalu- 
ated at the set of points (0, t) and (z, t'). In other words, 
we assume that the particle trajectory is not so strongly 
correlated with the surface slope configuration, such that 
averages of surface characteristics are not affected. Nu- 
merical results show that this approximation is justified 
except when the particle motion is slower than the sur- 
face fluctuations. We shall discuss this case shortly, and 
assume for the time being that this assumption is valid. 
This approximation is the most crucial step in our calcu- 
lation. 

For the Edwards- Wilkinson surface, the slope correla- 
tor can be computed exactly, and the result is (for t > t'), 



{d,h{o,t)d,h{z,t')) 



D 



v\v{t- t') 



-exp 



v{t-t'), 



(6) 



Let us now conjecture the following asymptotic scaling 
form for the probability distribution function of particle 
displacement: 



P(t,f',.) = i/(j 



(7) 



where ^{t—t') is the characteristic dynamic length scale 
of particle motion. The scaling function f{x) is required 
to satisfy the normalization f{x)dx= 1. It also fol- 
lows that the mean-square displacement is 



{x\t))=^,e 



(8) 



where /i = f{x)x^dx. We now substitute Eq.6 
and Eq.7 into Eq.5. The effective noise correlator thus 
becomes 



$(t,0 = ^ffl ■,T = \t-t'\ 



where 



G{T) 



f 

Jo 



vT 



(9) 



(10) 



We note that G{T) is proportional to the Laplace 
transform of the function g{x) = x~^ f{x), i.e., G{T) = 
^^ g(A) with A = We now conjecture a power-law 
growth of the characteristic length scale ^ at late times: 



m 



(11) 



where (f) is the dynamical exponent that characterize 
the particle motion and to is a microscopic time scale. 
From the preceding equations, it is clear that consistency 
requires (j) > \. For, if < ^, the l.h.s in Eq.3 would 
vanish at large t, whereas the r.h.s would not. The value 
= i is also inconsistent with this equation, since the 
l.h.s would be a constant whereas the r.h.s would grow 
as ^/t at large time t. Thus the only possible values are 
4> > ^. We now focus on the large-i limit, where the 
Laplace transform variable A ^ oo. The behavior of 
G{T) in this limit follows from the following theorem [8]: 

limx^ox^''g{x) ^ limx^oo ;p>-l (12) 

It is reasonable to assume that /(O) is a non- vanishing 
constant, in which case g{x) ~ as x 0. From 
Eq.l2, it then follows that p = — |, and hence ^(A) ~ 
as A — > oo. Thus, the asymptotic form of G{T) is 



2 



G(T) 



7rD/(0) 



oo 



(13) 



After substituting Eq.8 and Eq.9 into Eq.3, and using 
the large T-form in Eq.l3, the resulting self-consistent 
integral equation for ^{t) is 



?^=--/(o) 



(14) 



We now substitute the power-law scaling form Eq.ll 
for ^(T), which gives 



oo 



(15) 



After equating the powers of time on both sides, we 
find = |, which is the principal result of this paper, 

along with the pre-factor a ^ ("^^) ^ ■ ^® have thus 
arrived at the somewhat surprising result that the slid- 
ing particle moves faster than the surface fluctuations, 
whose correlation length scales with time as ^/id. This 
super-diffusive behavior formally has its origin in the 
long-range nature of the effective noise correlator $(t, t'), 
which, from Eq.9 and Eq.l3, has a power-law tail of the 
form \t — t'\^'^. Similar anomalous diffusive motion of ad- 
vected particles with </> = § has been shown to occur in 
the one-dimensional (non-linear) Burger's equation with 
noise [3,7], using a mean-field approach. 

Before proceeding to discuss the numerical results, we 
would like to point out that the calculation presented 
above is valid only if lirrit^ooj^ vanishes, else we would 
not be able to use Eq.l2 to arrive at our result. Indeed, it 
is quite possible that in the asymptotic limit, ^{t) ~ Vui, 
but that this regime lies outside the validity of this ap- 
proximation. In such a case, the anomalous diffusion at 
early times can be still shown to exist, but with a dif- 
ferent exponent. To see this, let us consider Eq.lO again 
at sufficiently early times t so that £,^{t) <C V^- In this 
case, G{T) ~ ^ is a constant, and after substituting in 
Eq.9, and subsequently in Eq.3, we find that £_{t) ^ t^^'^ 
when ^ <C \fvi. From numerical simulations (to be pre- 
sented in the subsequent section), we are unable to rule 
out this possibility. Further discussions on this point are 
presented at the end of Sec. III. 



III. NUMERICAL RESULTS 

In this section, we discuss the results of numerical sim- 
ulations of the problem. The surface is constructed as a 
set of height variables {/li}, where i = l,2,..iV and N 
is the lattice size. Periodic boundary conditions are im- 
posed on the surface. We use the single-step model of the 
surface [9], where the height difference between adjacent 
lattice points is restricted, i.e., = ±[hi±i — hi] = ±1, 
where we have defined the clockwise/anti-clockwise slope 



variables t~ . The surface evolves through the following 

set of dynamical rules: 

• If = 1 and T~ 

bility q{q< 1). 



-1, hi hi + 2 with proba- 



If r+ = 
bility q. 



-1 and = 1, hi ^ hi — 2 with proba- 



The large-distance, long-time properties of this model 
are identical to the continuum model in Eq.l. We start 
from a 'ordered height configuration, where /i^ = 1 is 
i is odd, and hi = when i is even. A lattice site is 
selected at random, and its height is updated in accor- 
dance with the above rules. For the surface configuration, 
one Monte Carlo step is counted after N such attempted 
updates, where N is the lattice size. The surface then 
evolves through T ~ lOA^^ Monte Carlo time steps in or- 
der to reach the steady state. We used two lattice sizes, 
N = 1024 and N = 4096 in our simulations. 

For a certain surface configuration in the steady state, 
we place 100 particles on the surface at randomly selected 
locations. At each Mcmte Carlo step, after a surface up- 
date is complete, a particle is selected at random (with 
initial position at lattice site i, say) and its position is up- 
dated according to the following set of dynamical rules. 

• If = 1 and = 1, the particle moves one lat- 
tice step in the clockwise direction with probability 
p{p< !)• 

• If = — 1 and T-~ = —1, the particle moves one 
lattice step in the anti-clockwise direction with 
probability p. 

Essentially, the rules say that the particle only slides 
down when the slope is favorable, and remains stationary 
when it is at a local minimum or maximum. The particle 
displacements are measured such that for every move in 
clockwise direction, X{t) ^ X{t) + 1, and for every move 
in anti-clockwise direction, X{t) ^ X{t) — 1. The parti- 
cles are assumed non-interacting, which is consistent with 
the single-particle picture. When the probability p > q, 
the particles moves fast relative to the surface motion and 
vice-versa. Also, without any loss of information, we set 
max(]3, q) = 1. The parameter we tune in the simulations 
is the ratio r = ^. The simulations were run up to 10'' 

Monte Carlo steps for N = 1024 and IQ^ Monte Carlo 
steps for N = 4096. We computed the probability dis- 
tribution of the displacement of the particles, as well as 
the mean square displacement, as functions of time. The 
results were averaged over 100 different starting config- 
urations for A'^ = 4096 and 1000 starting configurations 
for N = 1024. 

In Fig.l, we have shown the results for the RMS dis- 
placement of the sliding particle plotted against the quan- 
tity qt, where t is the number of Monte Carlo steps. The 
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fraction q is used to scale the time so that all the surface 
configurations pass through an almost equal number of 
updates (In the continuum language, this is equivalent 
to using the combination vt). We have shown the results 
for four values of the update frequency ratio r: 0.01, 0.1, 
1 and 10. The results for r = 0.01 and r = 0.1 offers 
evidence for the power-law behavior ^ ~ is as predicted 
by the self-consistent argument. For r > 1, on the other 
hand, the growth of the RMS displacement with time is 
much slower. In this case, the observed growth exponent 
is closer to 5 at all times. The case r = 1 displays in- 
termediate behavior, with an effective exponent close to 
0.56. 
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FIG. 1. The figure shows the time evolution of the 
root-mean-square displacement of the sliding particle for sev- 
eral values of the update frequency: r = 10, 1.0.1 and 0.01 
(top to bottom). The short straight lines are guides to the 
eye, and have slopes 0.5, 0.56, 0.68 and 0.67 (top to bot- 
tom). When r < 1, the dynamic exponent is close to 2/3, 
whereas when r > 1, it is close to 1/2. For r = 0.01, there 
is an carly-timo regime whore the growth exponent is smaller 
than 2/3, but the asymptotic behavior is identical to that of 
r = 0.1. The lattice size in these simulations is iV = 4096. 



Although it is tempting to say that the cases r > 1 
and r < 1 are characterized by different dynamical expo- 
nents, we would like to be cautious here. Our numerical 
data does not rule out the possibility that the observed 
anomalous diffusion of the sliding particle when r < 1 is 
simply a transient regime, and the asymptotic behavior 
might be identical for all values of r. In the scenario dis- 
cussed in the last paragraph of the previous section, the 
observed behavior might reflect a slow cross-over from 
^ ~ t^/^ to ^ ~ t2. The data for r = 0.1 shows some 
bending for t > 10^, but at the moment we cannot con- 
clusively say whether this shows the cross-over towards 
a diffusive regime, or simply a finite size effect. We have 
checked that a smaller system size (A^ = 1024) shows 
bending as early as i ~ 10^ (Fig. 2), which makes it likely 
that this is a finite size effect. 



FIG. 2. The figure illustrates finite size effects in the time 
evolution of the root-mcan-square displacement of the sliding 
particle for r = 0.1 (top) and r — 0.01 (bottom), for lattice 
sizes N = 1024 and N = 4096. The short straight Une at the 
top is a fit, and has slope 0.67. 

Particle trapping in the valleys: The origin of the 
distinct scaling behavior when r > 1 may be under- 
stood through the following argument. When the particle 
moves faster than the underlying surface fluctuations, af- 
ter the initial downward slide towards the local minimum, 
it gets 'trapped' there. At all future times, the motion of 
the particle is 'slaved to the dynamics of this local min- 
imum. The long-time dynamics of the particle may be 
visualized using the concept of valleys in the surface. A 
valley is defined as a linear stretch of the surface below 
a certain reference height, which we choose to be the in- 
stantaneous height at the starting location of the particle 
(Fig. 3). It is well-known that any instantaneous height 
configuration of a EW surface relative to a reference point 
on the surface is a random walk (RW), which returns to 
the origin after N steps. Thus, the length distribution of 
the valleys is simply the distribution of return times of a 
one-dimensional RW to its starting point, which follows 
a power-law decay: P{1) ^ l~2 for Z <C A" [10]. Let us 
now consider a section of the surface with length ^, and 
let N^{1) be the number of valleys of length I within this 
section. Since we require / dlN^{l)l = ^, it follows that 
N^{1) has the scaling form 

mi) ~ V^r^ (16) 

The slope correlations of the surface over a length I will 
decay by a time t{1) ~ I'^/v from Eq.6. For a valley of 
length /, this time is roughly the time scale for the over- 
turning of the valley, whereby the valley is transformed 
into a hill. An event of over-turning leads to a displace- 
ment of a particle trapped inside the valley by a distance 
~ I. The total time required to produce a displacement 
~ ^ is then the sum of the trapping times inside all the 
valleys within ^, which is given by 
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tiO ^ [ airmail) (17) 
Jo 

After using Eq.l6, this equation gives t{^) ~ C^/z^- 
Equivalently, for a fixed time t, the typical displacement 
of the sliding particle is ^(t) ^ ^/id. The result shows 
that in the regime r > 1, the dynamical exponent for the 
sliding particle is |. 




FIG. 3. Illustration of a valley in the single step model: The 
open circle shows the position of the particle at time t = 
and the filled circle shows its current position. The horizontal 
straight line divides the surface into valleys and hills. 




FIG. 4. The figure shows the scaled probability distribu- 
tion f{x) = {X^ (t))^ P{0,t, z) plotted against the scaled dis- 
tance X = rf^ {X'^if)) for three widely separated values of 
time. The good scaling collapse justifies the scaling form used 
in Eq.7. The lattice size is Af = 1024 and the update fre- 
quency ratio is r = 0.1 in this figure. The parameter a « 1.45 
in the Gaussian fit, and is found to be the same for all val- 
ues of r. The data represents an average over 1000 different 
surface starting configurations. 

Finally, we discuss the results for the probability dis- 
tribution P(0, i, z) of the displacement of the particle at 
time t, relative to its location at time t = 0. In Fig. 4, we 
have displayed the results for r = 0.1 at three widely sep- 
arated instants of time, t — 10^, 10"^ and 10''. When the 

scaling function j(x) = {X'^{t))'^ P{0,t, z) of the proba- 
bility distribution is plotted against the scaled distance 



X = zj {X'^{t)), we find good collapse of the data, which 
provides a a posteriori justification for using this scaling 
form. It is also remarkable that the scaling function is 
very well represented by a Gaussian (see the fit in the 
figure). 



IV. CONCLUSIONS 

To conclude, we have shown by means of a self- 
consistent argument that the motion of a passive sliding 
particle on a fluctuating surface may show non-trivial 
behavior even when the surface is Gaussian, like the 
Edwards- Wilkinson surface. Specifically, the effective dy- 
namical exponent for the motion of the sliding particle 
may be the same as, or different from the dynamical 
exponent of the surface fluctuations, depending on the 
relative time scale of surface and particle motion. The 
observed apparent super-diffusive dynamics of particle 
displacement in the EW surface urges us to be cautious 
when the measurements of the passive sliding particle is 
used to measure the dynamical exponent or other dynam- 
ical properties of the surface itself, as has been suggested 
in [3] . It might be interesting to extend this study to look 
into the effects of changing the relative update frequency 
of the particle and surface conflgurations on the steady 
state characteristics of the problem in a many-particle 
context [2,4]. Lastly, simulations on a larger time scale 
using much bigger lattices would be necessary to deter- 
mine conclusively if the observed anomalous diffusion is 
indeed a cross-over effect. 
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